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STUDY OF COAXIAL HALL CURRENT ACCELERATORS 
AT MODERATE PRESSURES 

By H. A. Hassan 
North Carolina State University 

Robert V. Hess and William Grossmann, Jr. 

Langley Research Center 

SUMMARY 

A study of coaxial Hall current accelerators at moderate pressures (of the order of 
10 mm Hg) and moderate mass flow rates (of the order of 1 g/sec) is presented. The 
analysis is based on the conservation equations of the plasma constituents and employs a 
Saha equation with the gas temperature replaced by the electron temperature. The arc 
voltage is shown to be a linear function of the applied magnetic field strength at constant 
total current; also, a voltage plateau is reached at the higher currents for a constant mag- 
netic field strength. In addition, for a given total current, magnetic field strength, pres- 
sure, and mass flow rate, the theory predicts the electron, ion, and neutral temperatures; 
degree of ionization; Hall current; and axial and azimuthal mean velocities. Comparison 
of theory with experiment shows that the theoretical model is in fair agreement with exper- 
iment. The influence of ionization rates and three-dimensional effects is discussed. 

INTRODUCTION 

Experiments with coaxial Hall current accelerators over a wide range of operating 
conditions (refs. 1 to 4) show that these devices have certain characteristics in common 
with the homopolar device (ref. 5). The voltage remains approximately constant, 
increasing current for a given magnetic field establishing a so-called "voltage plateau," 
and the increase in voltage with magnetic field for a given current is approximately linear. 
Since the ionization controls the voltage-current and the voltage — magnetic -field charac- 
teristics, the common characteristics could be a result of the fact that the field geometry 
in the ionization region of the coaxial Hall device is similar to that in the homopolar device. 

Even for this similarity in configuration, the theories proposed for the homopolar 
device, where the mass flow rate is zero, cannot be applied to the coaxial Hall current 
accelerator, as has been recognized in reference 6. The theory proposed by Alfven 
(ref. 5) and extended by Petscheck (ref. 5) and Lin (ref. 7) for the homopolar device is 



restricted to the case where the ions as well as the electrons orbit with azimuthal E/B 
drift velocity where E is the electric field strength and B is the magnetic field 
strength; thus, azimuthal Hall currents are not permitted to flow. (The possible azi- 
muthal current due to differences of centrifugal effects on ions and electrons is neglected.) 
It is assumed by Lin that under such conditions, the ion temperature is in excess of the 
electron temperature and the electrons gain their energy for ionization from the ions. 

This condition is not in agreement with temperature measurements in coaxial Hall current 
accelerators at moderate mass flow (ref. 1) as well as at low mass flows (ref. 2) where 
the electron temperature is considerably in excess of the ion temperature t^. 
Dobryshevskii (ref. 8) shows that excess of ion temperature occurs when u>. r. becomes 

equal to or greater than unity or in the case where the degree of ionization is large. 

(u>i is the ion cyclotron frequency and r- is the mean time between collision of ions 
and neutrals.) He bases his theory for the homopolar device on an extension of the theory 
of the positive column to the case of a discharge in a transverse magnetic field with 
o> e 7^ e » 1 and < 1 (o) e is the electron cyclotron frequency and is the 

mean time between collisions of ions and electrons), whereby the electrons gain their 
energy directly from the electric field. 

Efforts to explain the approximate voltage plateaus and the linear voltage magnetic 
field characteristics in coaxial Hall accelerators have been made previously by the pres- 
ent authors (refs. 1, 9, and 10) and by Patrick and Schneiderman (ref. 4). The approach 
used in reference 4 attempts to apply the theory proposed by Alfven and Petscheck for the 
homopolar device to the coaxial Hall current accelerator. In this theory the ion slip is 
assumed to be very large; that is, the interaction of ions and neutrals is negligible and the 
ions orbit with the E/B drift velocity. It will be shown in the section "Analysis" that 
the approach used in reference 4 must be modified if the ion slip is moderate and ions are 
put into rotation progressively along the axis. This modification should be especially 
important for the acceleration mechanism postulated in the conclusions of reference 4 
which are based on acceleration of the neutrals by a small fraction (5 percent) of ions. 
Furthermore, although the assumption is used that the ion slip is large for the analysis 
of the voltage -magnetic field strength characteristics, in reference 4 the analysis of 
the current voltage characteristics is based on considerations which are valid only if 
the ion slip is small. 

It is thus seen that theories for the homopolar device and that of reference 4 for the 
coaxial Hall device are based on assumptions that are not generally valid in a coaxial Hall 
current accelerator. The theory presented here (see also refs. 1, 9, and 10) does not 
depend on the assumption that the ions are orbiting with the drift velocity. It is based on 
the conservation equations of the plasma constituents that may be derived (for example, 
ref. 11) from the appropriate Boltzmann equations. The flow is assumed to be steady 
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and one -dimensional. This assumption implies that the effects of the electrode sheaths 
are small. Measurements of comparatively small power losses to the electrodes in 
reference 1 for moderately high mass flows and pressures support this assumption. The 
conditions under which the electron temperature is expected to be higher than the ion 
temperature are discussed in detail. It is concluded that for typical operations of coaxial 
Hall current devices, the electron temperature is always greater than the ion temperature. 
To keep the numerical complications at a reasonable level, the assumption of local ioni- 
zation equilibrium has been employed; Saha’s equation, the gas temperature being replaced 
by the electron temperature as suggested by Kerrebrock (ref. 12), is used. Actually, the 
residence time of particles for various operating conditions of coaxial accelerators is 
generally less than the time required to reach ionization equilibrium. Although appre- 
ciable deviations from ionization equilibrium still exist, even for the higher mass flow 
rates (ref. 1), the simplified theory still provides a reasonable agreement with experiment. 
Preliminary results using the rate equation were presented in reference 10. 

The calculations were carried out for two and three temperature models and various 
effective cathode radii. The dependence of the various collision cross sections on tem- 
perature was accounted for by using the best available experimental data. As a result of 
these calculations, the electron, ion, and neutral temperatures, the degree of ionization, 
the Hall current, the voltage -current, and the voltage-magnetic field characteristics can 
be predicted for a given magnetic field, total current (or power input), pressure and mass 
flow rate. It will be shown that regardless of the value of the effective cathode radius, 
the voltage was a linear function of the magnetic field at a constant current. Also, a 
voltage plateau is indicated for a constant magnetic field at the higher currents. The 
other predictions of the theory were in fair agreement with experiment. It is assumed in 
the analysis that the discharge is uniformly distributed. The transition from spokes to a 
uniformly distributed discharge has been studied in reference 1 as a function of magnetic 
field and temperature of the arc-heated thoriated tungsten cathode. 

The work done by H. A. Hassan presented herein was done under NASA Grant 
NsG-363 at North Carolina State University, Raleigh, North Carolina. The analysis pre- 
sented herein was completed in July 1965. 

SYMBOLS 


A cross-sectional area 

a cathode radius 

B magnetic field strength 

b anode radius 


3 


D 

E 

e 


Ss 

H 

h 

I 

i 

J 

i 

K 

k 

l 

m 

m 

n 

n s 

P 

P 

Q 

Qi 

q 

*s 

r,e,z 


operator 

electric field strength 
electronic charge 

degeneracy of the ground state of species s 

specific enthalpy 

Planck’s constant 

ionization potential 

unit vector 

total current 
current density 

quantity defined by equation (6) 

Boltzmann constant 

length of discharge 

mass of particle 

mass flow rate 

number density of nuclei 

number density of species s 

power input 

pressure 

rate of heat loss 

ionization cross section 

heat flux 

rate of production of species s 
cylindrical coordinates 
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T temperature 

t time 

tj characteristic ionization time 

t r characteristic residence time 

U mean velocity 

u,v velocity components in directions of z and 6 

V voltage 

V s velocity of species s 

Z s j. diffusion cross section 

a degree of ionization 

e energy of internal modes 

0 dimensionless temperature of species s 

s 

X quantity defined by equation (18) 

p density 

p_ density of species s 

s 

ct q conductivity in absence of a magnetic field 

r g £ mean time between collisions of particles s and 

o» s cyclotron frequency of species s 

Subscripts: 

i,e,a ion, electron, neutral 

s,t refer to any of the plasma constituents 

r,0,z refer to cylindrical coordinate system (fig. 1) 

An arrow over a symbol denotes a vector quantity. 



ANALYSIS 


The analysis is based on the conservation equations for the plasma constituents. 

In spite of the fact that the flow field in a coaxial Hall current accelerator is axisymmet- 
ric, with both r and z dependence, it will be assumed, as a first approximation, to be 
a function of z with r appearing as a parameter. For a steady, one -dimensional flow, 
the overall continuity equation can be written as 

puA = m (1) 

where p, u, A, and m denote, respectively, the mean density, mean axial velocity, 
cross-sectional area, and mass flow rate. The individual continuity equations for the 
plasma constituents may be written as 



where n g , V g , and R g are the number density, velocity, and rate of production in par- 
ticles per volume of species s, respectively. To avoid the unnecessary complications of 
dissociation and multiple ionization, only singly ionized monatomic gases are considered. 
For such gases, ionization is due mainly to electron-neutral collisions; that is, the domi- 
nant reaction is 


A + e — A + + e + e (3) 

Opposing this reaction is the recombination reaction 

A + + e + e — A + e (4) 


In equations (3) and (4), A denotes the atom; A + , the ion; and e, the electron. Based 
on these reactions, R g can be expressed as 


R - R e - Rj - -R a - 1 n e n a ” 


mn, 


l e 
K 


f(T e ) 


(5) 


where 



( 6 ) 


and (from ref. 13) 


- (4;) 1/2 (s;) 3/2 If “*(- 


3/2 


(7) 
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where gj,g a are the degeneracies of the ground states of the ions and atoms; m e , the 
mass of the electron; £, dummy variable, k, the Boltzmann constant; h, Planck's con- 
stant; I, the first ionization potential; T e , the electron temperature; and Qj is the 
ionization cross section. If the diffusion velocity is small compared with the mean veloc- 
ity, that is, if 

V g * U (8) 

where U is the mean velocity, equation (2) gives for a one-dimensional flow 

nug = R s (9) 

where a is the degree of ionization defined by 


n e = n. = an 


n a = (1 - a)n 


( 10 ) 


and n is the number density of the nuclei and where the gradient of the diffusion velocity 
is assumed to be small. 

Because of the complexity of equation (9), the question arises whether a simpler 
equation can be employed to calculate the degree of ionization. Examination of equa- 
tions (5) and (9) shows that a characteristic ionization time can be defined as 


<U) 

Thus, t^ is small if f(T e ) is large. This characteristic time may be compared with a 
typical residence time t r . For the case where 


h 

t r 


< 


1 


ionization equilibrium is attained in the chamber; in this case 


n n 
e a 


_ n i n e 


K 


or 


a 2 _ K 
1 - a n 


( 12 ) 


(13) 
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which is Saha's equation. On the other hand, if — > 1, ionization equilibrium is not 

tr 

attained in the chamber and equation (9) has to be employed. It is expected that for the 
higher mass flow rates and the higher pressures, equation (13) can be employed. This 
procedure was followed for the case of moderate mass flows for the sake of simplicity. 

The momentum equations give, in addition to the expressions for the mean veloc- 
ities, the ion and electron current densities. If the coupling between diffusion and heat 
flux is ignored, one can use the expressions derived by Cowling (ref. 14) for the ion cur- 
rent density and the total current density. Assuming that collision cross sections of 
charged particles are much larger than cross sections of collisions involving neutral 
atoms and solving Cowling’s expressions explicitly for the current densities yields 



(14) 

(15) 

(16) 

(17) 

(18) 




(19) 


( 20 ) 
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where E, B, and Z s t denote, respectively, the electric fields, magnetic field, and dif- 
fusion cross section. 

If viscous stresses are neglected, the overall momentum equation can be expressed 
as 


p (^) +vp= ^ ><5 ] 

P = Y Ps = Y kn S T s 


( 21 ) 


where p is the pressure. In the arc region of a coaxial accelerator, a crossed-field 
geometry with an axial magnetic field and a radial electric field is assumed. Thus, in the 
region under consideration, 


u = i z u+ i e v 


j = - i r j r " 1 03e 


( 22 ) 

(23) 


B= i z B 


E = -i r E 


(24) 


The components of equations (21) in the z- and ^-directions can be expressed as 


puJU^Uo 

dz dz 


(25) 


and 




(26) 


Also, use of equations (22) to (24) in equation (14) yields 


V 


j r = or(E - vB) 

a, e T ie( 1 + x2 ) 


(27) 


and 


i - 3r 

~T 


( 28 ) 
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If one chooses 


j r = 


J 

2jrr l 


(29) 


where J is the total current and l is the length of the region over which the discharge 
extends , the charge conservation equation 


V • j = 0 


(30) 


is satisfied. The independence of j r with respect to z was chosen to keep the analysis 
simple. If one had assumed instead independence of E with respect to z, the variation 
of j r with z would have to be considered. The total current J applied over a length 
l would then have to be expressed in the more general form than equation (29); that is, 

J = 2jrr r j r dz (30a) 

;0 

The limitations of the assumptions regarding the field geometry and the one- 
dimensional approximation can be seen from equations (26), (27), and (29). Equation (29) 
shows that the properties are functions of r in addition to being functions of z. This 
condition would indicate that for axisymmetric devices, the one-dimensional approxima- 
tion may be rather crude, especially in the regions where derivatives with respect to r 
cannot be ignored. Also, for constant cross-sectional area, equation (26) can be readily 
integrated to yield 


v = 



dz 

J r 


(31) 


where for the case of j r independent of z, using equation (29) yields 

v = - 

m 2ttt i 


(31a) 


It is indicated from equations (27) and (31a) that if constant current density j r is 
to be maintained with increasing z, the electric field would have to increase with z to 


E r with z is small in the pres- 


E 

prevent v from exceeding — . Since the variation of 

B 

ent ejq)eriments, the simplified analysis using j independent of z is restricted to 

E ^ 

values of v small compared with — . Experimental evidence obtained in reference 1 

B 

indicates that this approximation is good for the moderate mass flows. Since for low mass 
flow used in electric propulsion a wide variety of electrode configurations exists, which 
have different distributions of E r and j r , a variety of distributions is being studied in 
order to permit comparison with experiments. 
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Three energy equations are needed to determine the electron, ion, and neutral tem- 
peratures. These energy equations may be taken as the overall energy equation and two 
of the species energy equations. If viscous dissipation is ignored, the overall energy 
equation may be written as 

P^(h + ^ + V.q-E.J = 0 (32) 


where 


H - $1 Vs - ta e) 

is the specific enthalpy and q is the heat flux vector. The species energy equations 
can be expressed as (ref. 11) 


d_ 

dz 



s 



+ 






K st 


k(T s - T t ) 
m s + m t 


(33) 


where 



m s m t ... 
m s + m t | 


2k (^ + ^) Z st n s n t 


(34) 


where 


e a = 0 


e i = I 


e e, = I 


(e e =°) > (35) 

(e.=0) 


and q s is the heat flux vector of species s. Since the mass of the atom is almost equal 
to the mass of the ion, the assumption is usually made that the ion and neutral tempera- 
tures are equal. As a result of this assumption, only equation (32) and one of equations 
(33) are needed. 


To complete the formulation of the problem, expressions for the various collision 
cross sections are needed. The ion-electron cross section is obtained by assuming that 
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the ion-electron collisions follow Coulomb's law. In this case, however, the resulting 
expression depends on the cut-off distance. If the suggestion of Spitzer (ref. 15) is used, 
for argon, 


Z ie = 2.96 X l(T l6 0 e - 2 log t 


e. 


3.441 xlO 6 -£-] 
V n e 


3\l/2 


(36) 


where 



(37) 


The electron-neutral collision cross section for argon was obtained from the data 
of Massey and Burhop, Barbiere, and Frost and Phelps (refs. 16 to 18). For electron 
temperatures of the order of 1.0 eV, the cross section has a minimum as a result of the 
Ramsauer-Towsend effect; whereas for electron temperatures of the order of 10.0 eV, the 
cross section has a maximum. By using the data of references 16 to 18, the cross sec- 
tion was approximated by the quadratic expressions, 

Z ea = °- 16 x 10 " 16 + 5 ’ 015 x lo_14 ( 0 e “ O-O 1 ^) 2 (38a) 

and 

Z ea = 15.835 x 10" 16 - 32.390 X lO“ l6 (0 e - 0.812440) 2 (38b) 


The data on ion-neutral collision cross sections for argon have been summarized by 
Fay (ref. 19). For temperatures less than 1000° K, the polarization force is important; 
whereas for higher temperatures, short-range forces are dominant. Since it is anticipated 
that ion temperatures higher than 1000° K exist in discharges of interest, it is assumed 
here that ion-neutral collisions follow the Lennard-Jones potential. At the higher temper- 
atures, the repulsive part of the potential is dominant; therefore, for argon (ref. 20), 


,Je i+ flaY 1 / 6 

Z ia = 2.528 X 10 -15 l — 1 ) 


(39) 


It should be noted that values of Zj a based on charge exchange are more appro- 
priate. However, the results were practically unchanged for values one order of magni- 
tude higher than those in equation (39). 

In problems of arc discharges, it is desired to calculate the arc properties for a 
given power input or a given total current. Since 
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P = JV 


(40) 


1' 

1 


where V is the voltage, an expression for V is needed. Equations (27), (30), and (31) 
give 

E = VieU^J_ + AJB?z (41) 

°o X 2nrl m 2ttt l 

Since E z is assumed to be zero, V is given by 


V = 


fb 

\ E dr 
J a 


(42) 


where a is the effective cathode radius and b is the anode radius. Equations (41) and 
(42) show that V depends on z; therefore, the assumption that E z = 0 is an approxi- 
mation if J cannot vary with z. Also, since the flow properties depend on r implic- 
itly, integration of equation (42) requires the solution of the one-dimensional equations 
for different values of r. In view of the assumptions introduced here, such a procedure 
is not warranted. Therefore, it is assumed that the dependence of E on r is one of 
proportionality to hence, 


V 


'Vie 1 + X 2 J 
l °o X 2i tI 


A JB 2 \ 
m 2 t r / 



(43) 


In the evaluation of expression (43), z was chosen to be equal to l . 

Examination of equation (43) shows that the voltage depends explicitly on the three 
parameters a, b, and l. Of these parameters, only b is probably well defined; it is 
usually chosen as the actual geometric radius of the anode. Determining the effective 
cathode radius is a very difficult theoretical and experimental problem, and current 
bulging makes it very difficult to determine l . The situation is made more difficult 
because a and l vary with J and B. Since a one -dimensional theory cannot predict 
a and l , these parameters are to be treated as experimentally determined parameters. 

From equations (30) and (42), the volume integral of equation (32) can be written as 

= JV- - Q (44) 

where the subscript o denotes conditions at the chamber entrance and Q is the rate of 
heat loss. In obtaining the volume integral, the assumption was made that E r and V 
are independent of z; the total current J, as indicated in equation (30a), is, however, 
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independent of the variation of j r with z. Equation (44) is valid regardless of any 
other assumptions regarding ionization rates or any other assumptions employed in 
deriving equation (14). 

THE CONSTANT -PRESSURE APPROXIMATION 

The geometry of the channel and of the field configuration chosen here limits the 
Mach number to a maximum of unity. It is expected that in the regions where the Mach 
number is less than one, most of the energy supplied to the arc goes into ionization. 
Because a small fraction is going into acceleration, the assumption of constant pressure 
in the arc region is somewhat justified. As a result of this assumption, the axial velocity 
is calculated from equation (1) and is 


u = — = Constant (45) 

PA 

The electron energy equation can be simplified from the consideration that the elec- 
trons lose very little of their energy in a given collision. This condition would imply 
that terms containing the electron temperature gradient may be considered to be small 
compared with the power input and collision terms. Introducing this assumption into 
equation (33) and assuming negligible effects of inelastic collisions yields 

- - v k(T e - T t ) 

E * • le = »2/ K et m (46) 

The Two -Temperature Analysis 

Because the mass of the atom is almost equal to that of the ion, the assumption is 
usually made that the ion temperature is equal to the atom temperature. As a result of 
this assumption, two energy equations are needed and these equations may be chosen as 
equations (44) and (46). Choice of equation (46) as the second energy equation shows that 
the electron temperature is always greater than the atom or neutral temperatures. This 
condition is not in agreement with the ionization mechanism advanced by Lin (ref. 7), in 
which he assumed an ion temperature in excess of the electron temperature. 

It is seen from equation (33) that in the situations where pressure and temperature 
gradients and inertia terms are small compared with power inputs and collision terms, 
choice of the ion energy equation, which in this case reduces to 
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would lead to an ion temperature in excess of the electron temperature, a result that is 
in agreement with Lin's assumption. As a result, the assumptions T- = T and steady 
uniform (gradients negligible) plasma are not compatible. Also, the choice of equation 
(46) and not equation (47) as the second energy equation should be justified. 

The choice of the appropriate energy equation can be inferred from a comparison of 
the terms containing the temperature gradients and the power input terms in equation (33). 
As has been pointed out, the electrons lose very little of their energy in their collisions 
with the ions and atoms. On the other hand, because the mass of the ion is almost equal 
to that of the atoms, the ions lose an appreciable part of their energy when they collide 
with the atoms. Thus, since collision of electrons with the heavier particles is charac- 
terized by longer relaxation times, it is expected that |DT e /Dt| < |DTi/Dt|. The power 
received by the ions and electrons can be calculated by employing equations (14) and (15). 
For the configuration under consideration, 

7i • E* = ^ S - (a >e T ie )- 1 ]E*2 

w e T ie( 1 + x > 

and 

j e • E* = — E* 2 (48) 

("eV 2 * 1 + x2) 

Hence, 

J. • E* 

J -=r = 2(1 - a) 2 o: e T ie Wi T ia (49) 

3 e • E * 

Thus, in situations where the ratio (eq. (49)) is of order unity or less, equation (46) is the 
appropriate equation. On the other hand, if the ratio (eq. (49)) is much greater than unity, 
use of equation (46) may be in error. In such case, use of equation (47) is more appro- 
priate provided that terms containing dTj/dz are small compared with jj • E*. 

Most operating conditions of existing coaxial Hall current accelerators are such that 
2(1 - of)^o)^T^ a co e Tj e < 1. Hence, use of equation (46) is justified in this work. Since this 
discussion did not take into consideration inelastic collisions and radiation losses, these 
conclusions are expected to be valid for electron temperatures of the order of 2.0 to 3.0 eV 
or less. 

The Three-Temperature Analysis 

The relation between the ion and electron temperatures in a steady and uniform 
(negligible gradients) plasma can be investigated by utilizing equations (46) and (47). 

15 



( 50 ) 


Solving these equations for T e - T a and Tj - T a yields 


T e - T a = 


n i Z ie E * • j +6n a Z ia E* • 3, 
/3A 


and 


n Z- E* 
e le 


T + n 0 Z 0 Q E* 

a ae 
/3A 


(51) 


where 



and 


A = n i n a Z ie Z ae + 6n a Z ia( n i Z ie + n a Z ae> 


(52) 


It is seen from equations (50) and (51) that T e > Tj if 



or, if equations (49) and (52) are used, equation (53) can be written as 


l/m_\ 1/2 ^T+^ 

4\m e ) y T e 


> 2(1 - <x) 2 “iTia“e T ie 

Zae 


(54) 


since Zi a /Z ae ~ 10 2 and T e > Tj for all practical purposes. This analysis gives 
added support to the considerations justifying the use of equation (46) in the two- 
temperature analysis. This result is limited to the case where WjT ia £ 1. It should 
be noted that this result is independent of the magnitude of v. 
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RESULTS AND DISCUSSION 


The plasma properties at a given section can be determined from the simultaneous 
solution of equations (13), (44), and (46) in conjunction with equations (31), (43), and (45). 
For the three-temperature analysis, these equations have to be supplemented by equa- 
tion (47). The results of these calculations for a given total current J, magnetic flux 
density B, mass flow rate and pressure give the electron, ion, and neutral temperatures, 
the axial and azimuthal velocity, the voltage, and Hall current. 


As has been pointed out, the quantities r, a, b, and l are constants which have 
to be specified for any given calculation. In this work, various values of r and a have 
been chosen; however, no effort has been made to vary b and l. Since most of the 
ionization is expected to take place near the plane of the cathode tip, a does not corre- 
spond to the radius of the cathode cylinder. (See fig. 1.) The values of a selected 
ranged from constants independent of J and B to parameters depending on B. In 
particular, it has been found that, for this investigation, a choice of a as 


a = 0.02 



cm 


(55) 


over the range of B from 3000 gauss to 9000 gauss results in a reasonable fit of the 
experimentally measured voltage. The values of r investigated consisted of constants 
and values determined by the relation 


with a given by equation (55). By using various values of r, one can assess, in a very 
crude way, the importance of three-dimensional effects. 

The calculations were carried out for the accelerator shown in figure 1. The device 
has an anode radius b = 1.5204 cm and the cathode cylinder has a radius of 0.899 cm. 

In the calculations reported here, the mass flow rate was kept constant at 1.8 g/sec, 
the pressure varied from 10 to 19 mm Hg and l was chosen as 1 cm. All the calculated 
values shown in the figures correspond to Q = 0 and the station z = l where the meas- 
urements were made. Heat losses were ignored because the measurements of refer- 
ence 1 showed that they were small compared with the power input. 

The predictions of the three-temperature and the two-temperature analyses were 
almost identical. This similarity may be attributed to two reasons. As can be seen from 
figure 2, the calculated T i and T a are very close to each other. Also, the range of 
input parameters considered here was such that the calculated value of 
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2(1 - a!)^a> e Ti e co a Ti a was always less than 1. It may be recalled that this criterion was 
derived in trying to justify the use of equation (46) as the proper second energy equation 
in the two-temperature analysis. 

Figures 3 and 4 show a plot of the variation of voltage with magnetic field intensity 
for various values of a and r for a total current of 100 amperes. The figures show that 
regardless of the values of a or r, the voltage increases linearly with B at constant 
total current. The various interpretations for the linear V - B relationship in homo- 
polar devices and coaxial Hall current accelerators were reviewed in the introduction. 

It was concluded then that the assumptions employed by the various investigators were 
not applicable to the coaxial Hall current accelerator. It may be recalled that the analysis 
presented here assumes that the electrons receive most of the energy from the electric 
field with the result that the electron temperature is in excess of the ion temperature in 
accordance with experiment. Also, Hall currents, which play an important role in coaxial 
Hall current devices, are included. In an effort to find out the reasons for the results 
shown in figures 3 and 4, attention is directed to equation (46) and references 9 and 10. 
When equations (19), (20), (34), and (48) are used, equation (46) can be written as 



(57) 


The results of the calculations showed that for the range of J and B considered, 

100 amperes to 300 amperes and 3000 gauss to 10 500 gauss, respectively, v < E/B, 

6^ < 0 e , and X and 6 e are "weak" functions of B at constant J. Thus, approxi- 
mately, equation (57) can be written as 

| ■ f ki 1 + + x2 >»e - (c * Constant) (58) 


where, for the conditions analyzed in this work, c was less than unity. The fact that X 
is a weak function of B is a result of taking into consideration the ion slip term; other- 
wise, the validity of equation (58) would require o> e r ie » 1 or X « 1. Equation (57) 
shows also that the limit of E as B approaches zero exists. Denoting this limit by 
E 0 yields 
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ie 


3 1 + 


1 '“!“)(•. -v 


a 


j ei 


( 59 ) 


The fact that the linear 


V - B characteristic of constant current exists for 
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(c < 1) 


(60) 


v 



shows that the line ar V - B characteristic is not dependent on the ions drifting with the 
critical velocity ^2 ^ = |i as indicated in reference 4, based on discussion of the homo- 
polar device in references 5 and 7. 

At the lower mass flow rates, the quantity c could conceivably be of order unity 
because of the high electron temperature observed. This condition would explain the 
experimental results of the low mass flow devices observed in references 2 and 4. As 
stated at the end of the section "The Two-Temperature Analysis" for electron tempera- 
tures in excess of 2.0 eV to 3 eV, the effects of inelastic collisions and radiation losses 
would have to be included. 

The calculations for moderate mass flows have been also carried out with an elec- 
tron energy equation which took into consideration the heat loss by radiation and the 
inelastic collisions. Both of these contributions were found to have negligible effects on 
the results. The possible influence of losses to the electrodes is being evaluated in detail. 

Figure 3 also shows the sensitivity of the voltage to the effective cathode radius. 

In this connection, figure 3 suggests a method for estimating the effective cathode radius. 
The procedure consists of employing this analysis together with the experimentally deter- 
mined voltage and assuming a as an unknown to be determined from the simultaneous 
solution of the governing equations. Figure 4 shows that the r dependence of V is not 
as has been assumed in deriving equation (43). It also illustrates the importance of the 
three-dimensional effects in coaxial devices. 

The question arises (also discussed in ref. 8 for the homopolar device) concerning 
applicability of the theory to light gases where the ion cyclotron radius could be smaller 
than the electrode gap. For a given mass flow, however, the particle densities will be 
considerably higher for the lighter gas, a condition which reduces r^ a and somewhat 
compensates for the effect of increased uq on When ion-neutral collisions are 

rare, the electrons can still gain energy directly from the electric field over distances 
smaller than the ion cyclotron radius due to gradients in the potential distribution, whereby 
the electron conduction may become turbulent. Such a situation exists, for example, for 
the ionization process in the linear ion Hall accelerator and could also find application for 
certain range of operations in the present device. Careful distinction must be made 
between the role of given characteristic parameters for the ionization region and for the 
region where the acceleration occurs. Thus, especially for the lighter gases, the present 
ionization mechanism using T e > T^ may determine the E/B characteristics, and the 
condition of r Li smaller than the electrode gap may play a role in the acceleration. 

This problem has to be explored experimentally. 


19 


I 



Figure 5 shows a plot of voltage -current characteristics. At the higher currents, 
both theory and experiment show the existence of a voltage plateau. The theory shows a 
slightly positive characteristic at the lower currents when the experiment shows a 
slightly negative characteristic. An explanation for this discrepancy can be attributed to 
the use of Saha's equation (eq. (13)). Use of Saha's equation with the gas temperature 
replaced by the electron temperature gives an optimistic estimate of the degree of ioni- 
zation. Consideration of the various terms in equation (44) shows that, at a given section, 
it can be approximated by 


JV 

m 



or 


V 



(61) 


The faster rate of increase of a with J at constant B (see fig. 6) accounts for the 
slightly positive characteristic predicted by the calculations. It is expected that use of 
the rate equation, equation (9), in conjunction with the other conservation equations may 
correct this situation. In this case, the relaxation effects result in a slower increase of 
a with J at constant B. 

This effect may also be attributed to the fact that heat-transfer losses were ignored 
in the calculations. Such losses would tend to reduce the temperature and, in turn, would 
reduce a. The voltage-current characteristics must be consistent with those obtained 
from equation (43) based on the generalized Ohm's law and the momentum equation. The 
advantage of using the total energy equation is its more general character; however, equa- 
tion (43) shows dependency on certain parameters which are of great importance in deter- 
mining the operation of the device. 

The experimentally determined Hall current (fig. 7) shows that the Hall current has 

j 

a maximum. Since jg = JL, this behavior is expected because X consists of two terms, 
one proportional to B and the other to 1/B; X is minimum when 


w e T ie 


= 2(1 - a) 2 u iTi 


la 


(62) 


The calculations show that the maximum occurs at a higher magnetic field. This condi- 
tion may be due to the fact that the calculations underestimate the electron temperature 
(fig. 8). Higher electron temperatures increase co e Tj e and help bring about the equality 
in equation (62) at lower magnetic fields. 
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The possible influence of the measurement technique (see appendix) on the discrep- 
ancies between theory and experiment is now discussed. As the magnetic field increases, 
co e T i e becomes larger and the electrons are "tied" more strongly to the magnetic lines. 
Then the current pattern may bulge further downstream and, as a result, using the same 
cross-sectional area to calculate from the total Hall current may introduce some 
error. Also, if the total current would extend further downstream, the local Hall current 
probe would only register part of the total current, since the current regions removed 
from the probe give a smaller signal. Of course, other changes in current patterns also 
may occur with increasing magnetic field (see also, ref. 4), including the possibility of 
instabilities that would reduce the Hall current. 

Figure 6 shows that the calculations result in a degree of ionization slightly higher 
than the measured value. As has been explained earlier, use of Saha’s equation with the 
gas temperature replaced by the electron temperature (eq. (13)) leads to optimistic pre- 
diction of the degree of ionization. 

In spite of the good agreement between theory and experiment one cannot conclude 
that the use of Saha's equation to calculate the degree of ionization is fully justified. The 
calculated electron temperature is less than the experimentally measured value. The 
good agreement between theory and experiment is a result of using this lower temperature 
in an equation which predicts excessive values for the degree of ionization. The fact that 
basic trends like the voltage-current characteristics and the voltage-magnetic field char- 
acteristics are obtained in spite of these and other approximations suggests that they do 
not depend on the detailed mechanism of ionization. These results shed some light on the 
reason for the existence of similar trends of these characteristics for a great variety of 
operating conditions, as stated in the "Introduction." 

In spite of the large scatter of the experimental data (ref. 1), it can still be seen 
from figure 6 that the rate of increase of a with J is less than that predicted by the 
theory. This result lends support to the discussion regarding the voltage plateau. 

Figure 8 compares the calculated electron temperature with that measured by 
Langmuir probes for J = 100 amperes. The theory underestimates the electron temper- 
ature because ionization rates were ignored in the calculations. The interpretation of the 
measurements of reference 1 showed a discrepancy between probe and spectroscopic 
measurements. The spectroscopic measurements were made at a point inside the elec- 
trode region, whereas probe measurements were made near the tip of the cathode, outside 
the electrode region. Because of this, some discrepancy is expected. However, a recent 
reevaluation of the spectroscopic data of reference 1 at the Langley Research Center 
brought the spectroscopic measurements in line with the probe measurements. When 
operating at 200 amperes in the presence of argon lines AI and AH, the AH lines being 
more intense, an electron temperature in the range of 14 000° K to 20 000° K has been 
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estimated. This value is in good agreement with the probe measurements at 200 amperes, 
which ranged from 14 000° K to 18 000° K. Further discussion of the spectroscopic 
measurements is given in the appendix. 

CONCLUDING REMARKS 

The theory presented here is based on the conservation equations of the plasma 
constituents and does not assume that the ions are orbiting with E/B drift velocity 
where E is the electric field strength and B is the magnetic field strength. It pre- 
dicts the linear voltage -magnetic field strength characteristics at constant current 
observed by the various investigators and indicates the existence of a voltage plateau for 
a constant magnetic field at the higher currents. 

It is shown that an ionization mechanism based on electron temperatures in excess 
of ion temperatures is responsible for the observed characteristics. This result is in 
contrast to the mechanisms proposed by Alfven and used by Lin, Patrick, and 
Schneiderman. 

The model employed here predicts fairly well the gross characteristics of the 
coaxial Hall current accelerator. 

It is evident from the discussion, however, that relaxation effects play an impor- 
tant role in crossed field discharges. This conclusion is reached in spite of the fact that 
various effects seem to counteract each other with the result that reasonable agreement 
is obtained by using Saha's equation. It is concluded that, for a better understanding of 
the ionization and acceleration processes, a more elaborate analysis, which takes into 
consideration relaxation and three-dimensional effects is required. The fact that for the 
approximations made, the basic trends, like the voltage-current characteristics and the 
voltage -magnetic field strength characteristics, . were obtained in spite of these and other 
approximations is, however, in itself a significant result. The result suggests that these 
basic trends do not depend on the detailed mechanism of ionization and distribution of 
power in the plasma. This behavior helps to explain the existence of similar trends of 
these characteristics for a great variety of operating conditions. 

North Carolina State University, Raleigh, N. C. 
and Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 8, 1966. 
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APPENDIX 


DISCUSSION OF SPECTROSCOPIC AND 
HALL CURRENT MEASUREMENTS 

Spectroscopic Measurements 

Electron temperature determinations in the coaxial accelerator were made from 
identification of the emitted spectra lines rather than from relative line intensities. 

Since spatial nonuniformities exist in the plasma, and spatial resolution could not yet be 
achieved, complete knowledge of where the emitted spectral lines originate is not possi- 
ble. Thus the relative line intensity method may give an incorrect temperature based on 
lines from different regions. The line identification method, however, yields correct 
information of the range of temperatures over the entire radial region. For a spatially 
uniform plasma, both methods should produce the same results. 

Both methods of electron temperature measurement are valid only in a steady 
state. Calculations indicate that the residence time of a particle in the accelerator is 
shorter than the time required to reach steady-state ionization. Thus, if the equilibrium 
state can be expected to occur at higher T e , the electron temperature measurements are 
a lower limit, since the temperature is still increasing as a typical plasma volume ele- 
ment leaves the region under observation. 

Three radial positions of the electrode region were observed in the measurements. 
Access to the electrode region was achieved by drilling three parallel holes through the 
anode body and anode ring and inserting three stainless- steel tubes through the body up to 
the anode ring. With the line of sight through the tubes and the holes in the anode ring, 
the three radial positions between the anode and cathode could be observed. These radial 
positions were (1) a position just off the cathode surface, (2) one just off the anode wall 
surface, and (3) one which is half the distance between the cathode and anode walls. 

Spectrograms were obtained from each of the three regions at various magnetic 
fields and arc currents, and the films were scanned with a microdensitometer. The traces 
showed that the most intense lines were emitted from the center tube and the weakest lines 
from the region nearest the anode wall. The lines were identified and the electron tem- 
perature was determined from the presence of the lines by two methods. In one method 
the Saha equation, corrected for low density conditions, was used to find the variation of 
the relative population of the argon ground states with electron temperature. In the other 
method, the intensities of several observed excited AI and AH lines compared with 
electron temperature tabulated by Oertel and Raiford were used. Although the method 
using the excited lines is more accurate, reasonable agreement was found with the ground 
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state method, and the range of electron temperatures was found to be from 14 000° K to 
20 000° K, when operating at 200 amperes. As mentioned previously, this limit is the 
lower limit of temperature. The upper end of the temperature range appears to be in the 
center position with the temperature decreasing toward the anode wall and the cathode. 
Determination of the true radial temperature distribution would require use of an Abel 
inversion. 

For the electron density measurements, Stark broadening of the line was used. 
This measurement yielded only an order of magnitude estimate of the density. This 
method is suited to electron densities of 10 ^ to 10^ electrons per cubic centimeter, 
whereas at low densities, around 10^ electrons per cubic centimeter, the plasma line 
profile is of the order of the instrument profile, and thus interferes with accurate meas- 
urements. It is estimated that the accuracy of the electron density measurements used 
in this method are not much better than ±50 percent. 

Hall Current Measurements 

The Hall current was measured by placing a search coil around the plasma stream, 
the coil being embedded in a boron nitride nozzle immediately adjacent to the copper 
anode ring. The Hall current signal is obtained by "crowbarring" the arc current into a 
resistor load, and thereby obtaining a clean arc shutoff. The collapsing Hall magnetic 
field then induces an electromotive force and current in the coil; the charge obtained by 
integration of this current in a ballistic galvanometer is proportional to the Hall current 
in the plasma. The calibration of the Hall current was performed by placing a small 
many-turn coil in the approximate position where the Hall current was thought to be 
largest. 

The Hall current measurements were obtained in amperes, and the Hall current is 
assumed to be distributed in a single-turn loop. To permit comparison with the theoret- 
ical values of Hall current density j e the Hall currents were divided by a constant 
cross-sectional area with the height of the 0.76-cm electrode gap and with an assumed 
discharge length of 1 cm. 
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Figure 1.- Schematic diagram of accelerator. 
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Figure 3.- Variation of voltage with magnetic field for given r. J = 100 amperes; r = 1.2. 
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Figure 4.- Comparison of theory and experiment. Variation of voltage with magnetic field for varying r. J = 100 amperes. 
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Figure 7.- Comparison of theory and experiment. Variation of Hall current with magnetic field. 
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